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Abstract 

Electric field gradients and chemical shielding tensors of the stable mono- 
clinic crystal phase of ethanol are computed. The projector-augmented wave 
(PAW) and gauge-including projector- augmented wave (GIPAW) models in 
the periodic plane-wave density functional theory are used. The crystal data 
from X-ray measurements, as well as the structures where either all atomic, 
or only hydrogen atom positions are optimized in the density functional the- 
ory are analyzed. These structural models are also studied by including the 
semi-empirical Van der Waals correction to the density functional theory. 
Infrared spectra of these five crystal models are calculated. 

Keywords: Ethanol crystal, NMR, NQR, Density functional theory, 
GIPAW, PAW, Infrared spectra 



1. Introduction 

Organic crystals show many interesting phenomena and they are impor- 
tant in various applications The solid ethanol exists in several polymor- 
phic phases [2]. The stable crystal is monoclinic jsf. In addition, three glass 
phases exist: a plastic crystal, an orientational glass, and the true structural 
glass jil,^!- The glass phases have been investigated by X-ray diffraction jsj, 
Raman js] , dielectric relaxation , calorimetric measurements j^, 0] , neutron 
scattering 0], EPR jsj, and NMR relaxation measurements 

The ethanol crystal belongs to the Pc space group and melts at 159 K 
jsj. This crystal consists of zig-zag molecular chains which are connected 
by hydrogen bonding. Two types of molecules are present in each chain: 
one is with a C-C-O-H torsion angle of -63°, and other is a conformation 
with a torsion angle of 179°. It is not easy to grow good quality monoclinic 
crystals because of complex thermodynamic properties of the solid ethanol. 
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The ethanol molecule is known as one of the first systems for which NMR 
measurements in organic solvents were carried out [lo|. The temperature 
dependence of the relaxation time in various phases of ethanol, including the 
monoclinic solid, have been studied by NMR |9| . In spite of intensive nuclear 
resonance studies of various organic materials and a substantial improvement 
of experimental techniques, we have not been able to find any report of mea- 
surements for chemical shifts and quadrupolar parameters in the monoclinic 
solid phase of ethanol. 

Advances in magnetic resonance experiments have been followed by a 
development of the first principle calculations of corresponding parameters. 
These theoretical approaches help in analysis of often complicated NMR 
spectra by connected them to the underlying microscopic structure of mate- 
rials. Recently, the density functional theory (DFT) based gauge-including 
projector- augmented wave (GIPAW) method has been developed lll-ll4|. In 



a comparison with older quantum chemistry techniques, which use the clus- 
ter approximation, GIPAW is particularly suitable for computations of the 
magnetic response of crystals. This method works with periodic boundary 
conditions and therefore correctly describes extended periodic systems, for 
example lattice effects and non-covalent interactions often present in molecu- 
lar crystals. In addition, GIPAW is an extension of the projector augmented 



wave (PAW) method [l5|. PAW and GIPAW use special pseudopotentials 
which, by reconstructing all-electron density close to atomic nuclei, correctly 
account for electrons in this region. PAW method has been successfully ap- 
plied in calculations of nuclear quadrupolar parameters in the density func- 
tional approach 16|, [17 . 

It has been found that the GIPAW method produces adequate results for 
a wide range of organic solids 
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In this study we present a DFT 
plane wave calculation of the ethanol crystal. We calculate corresponding 
hydrogen, carbon, and oxygen chemical shieldings in the GIPAW model. 
Quadrupolar coupling constants and asymmetry parameters for and ^^O 
are computed using the PAW method. 

In contrast to NMR and NQR parameters, measurements of infrared spec- 
tra of the ethanol crystal have been carried out 2^, |23||. Infrared spectra of 



the ethanol crystal, as well as of thin and thick solid ethanol films deposited 
on several substrates, have also been measured because of applications in 
astrophysics 2J, |25|]. Using density functional perturbation theory methods 



26[, we calculate infrared spectra of the ethanol crystal and compare these 



results with experiments. 
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2. Computational methods 



NMR and NQR parameters of the monochnic crystal of ethanol are com- 
puted using the GIPAW and PAW models llNl4| implemented in the Quan- 
tum ESPRESSO suite of computer codes for electronic-structure calcula- 
tions and materials modeling at the nanoscale 27|. The crystal structure 



and atomic positions are taken from the Cambridge Structural Database 
28 1 where data measured by the X-ray diffraction at 87 K are available. 



Electronic ground state and NMR calculations using directly experimental 
atomic positions have been carried out in similar DFT studies for crystals 
whose structures are well known and used [l^. We take this approach here 
in our first crystal model (labeled as a in the results), but also studied sev- 
eral other geometries. In the second model (6), we take the same data from 
the Cambridge Structural Database, but then we optimize all atomic posi- 
tions in DFT within the fixed unit cell. A such approach is usually taken 
for older experimental data obtained by X-ray diffraction measurements in 
which positions of atoms were often determined with large errors. In the 
model c, when relaxing all atomic positions we added the semi-empirical Van 



der Waals correction at the DFT-D2 level [29|, |30[. This correction has been 
tested for large benchmark sets of molecules, clusters, surfaces, and crystals 
where the Van der Waals interaction is important, and good agreement with 
available experimental results have been found. It is known that the largest 
errors in older X-ray measurements of organic materials exist for H atom 



positions [12[. Therefore, we also study the models d, where only hydrogen 
atoms are relaxed in DFT, and e where only hydrogen atoms are relaxed in 
DFT corrected by the Van der Waals interaction. 

The generalized gradient approximation (GGA) of Perdew, Burke, and 



Ernzerhof j3l|, and norm-conserving GIPAW pseudopotentials prepared by 



A. P. Seitsonen [32| are applied. The cutoff energy of 90 Ry is used. The 



Brillouin zone is sampled with the 12x10x8 Monkhorst-Pack set of k points 



33|. We found that the f-sum rule is fulfilled under these conditions. As 
reference systems we study the isolated tetramethylsilane (TMS) and ethanol 
molecules positioned in a box of 40 Bohr. For these systems we also use the 
cutoff of 90 Ry." 

— * — » 

The external magnetic field B and an effective field Bin{r) in a material 
are connected by 

S,„(r) = -a(r-)5, (1) 
where a{r) is the nuclear shielding tensor. Chemical shifts are often calcu- 



3 



lated from 

5(r) = cr,,;-a(r), (2) 

where aref is defined for a chosen reference nucleus, and a{r) = l/3Tr[a{f)]. 
The parameters typically used to describe quadrupolar interactions are the 
coupling constant Cq and the asymmetry parameter 77. The quadrupolar 
coupling constant is 

Cq = (3) 

where e is the electron charge, Vzz is the z principal component of the EFG 
tensor, Q is the nuclear quadrupolar moment, and h is the Planck constant. 
The asymmetry parameter is 

Vq = 7^^, (4) 

''zz 

where the components of the EFG tensor satisfy the condition \ Vzz\ > \Vyy\ > 
\Vxx\- Quadrupolar moments for nuclei we study are: Qi^H) = 2.860 x 
10-31 m^ and Qi^'^O) = -25.58 x lO'^^ m^. Q- Using these values of 
nuclear quadrupolar moments and Eqs. (l)-(4) we calculate NMR and NQR 
parameters of the ethanol crystal in the GIPAW and PAW models. 

Infrared spectra of five models of the ethanol crystal are calculated to 
study the stability of lattices and to compare results with experiments. The 
methods based on a density functional perturbation theory [26| and the 
Quantum ESPRESSO package ji^] are applied. The same GIPAW pseu- 
dopotentials as in NMR/NQR calculations are used. We have been increas- 
ing thresholds for self-consistency in ground state and phonon calculations 
(to 10-1^ Ry and 10"^°, respectively), as well as the cutoff to 110 Ry, in 
order to check the ability of crystal models to get rid of negative frequencies 
of normal modes. 



3. Results and Discussion 

Table 1 shows bond lengths for different strategies in the optimization 
procedure. The average value is taken where more bonds of the same type 
exist. The 0-C and C-C bonds are with similar lengths in all models of 
atomic positions and agree at the 0.01 A level. Differences in bond lengths 
with H atoms are larger. The C-H, 0-H, and H bonds in methyl groups in- 
crease in DFT optimized models in a comparison to the experimental data. 



4 



The lengths of hydrogen bonds decrease in DFT-optimized models. It is 
possible to explain differences between experimental and optimized atomic 
positions by always present problems in DFT (for example, approximations 
in exchange functional and pseudopotentials), but also by the fact that X- 
ray measurements jsl are rather old. Therefore, as in other similar studies, 
positions of atoms are sometimes determined with large errors. In addition, 
experiments were carried out at 87K, whereas DFT calculations neglect ther- 
mal effects. Differences between DFT-optimized bond lengths without and 
with the van der Walls corrections (i.e., bvs c and dvs e sets) are nonexistent 
or small. The largest change is at the 0.01 A level when all atomic positions 
are optimized (6 vs c ), and at 0.001 A when only H atoms are relaxed {d vs 
e). However, it is known that calculations of NMR parameters are sensitive 
even to small changes of atomic positions. 

The chemical shieldings are shown in Table 2, and the corresponding 
atomic labeling is explained in Fig. 1. The results show that the largest 
shieldings exist for hydrogen atoms in the OH group, then for ones from the 
methyl group, and the smallest are for H in the CH2 group. Hydrogen atoms 
in the same methyl group have slightly different shieldings. Carbon atoms 
from the methyl group exhibit larger shieldings than C in the CH2 group. 
Shieldings in DFT-optimized models are much smaller then ones calculated 
for positions from the X-ray measurements. It is important to point out 
that our calculations do not consider temperature effects. NMR and NQR 
calculations are carried out for static molecules in the crystal lattice and for 
fixed atomic positions. 

The results of density functional theory calculations are the components 
of the nuclear shielding tensor. Different approaches are taken in the litera- 
ture to compare these results with experiments. Experimental results for the 
ethanol crystal are not available. Therefore, we are not able to fix the scale 
using the fitting to the experimental data. Instead, calculated here shield- 
ings in isolated molecules of TMS (for hydrogen and carbon), and ethanol 
(for oxygen) are taken as reference values: cr{H) = 30.4 ppm, cr{C) = 178.2 
ppm, cr(0) = 258.6 ppm. We did not use average values for a{C) and <j{H) 
calculated in the ethanol molecule because of the fact that these atoms exist 
(in a molecule, as well as in the crystal) in different geometrical and chemical 
environments. In contrast, shieldings are the same for all four carbon and 
all twelve hydrogen atoms in TMS. Because of the lack of the experimen- 
tal fixing scale, and to benefit possible future experiments, we report both 
shieldings (Table 2) and shifts (Table 3). 
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Table 1: Bond lengths (A) in the ethanol crystal, a: X-ray atomic positions b: all 
atomic positions are optimized in DFT, c: all atomic positions are optimized in DFT 
corrected by the van der Waals interaction, d: positions of hydrogen atoms are optimized 
in DFT, e: positions of hydrogen atoms are optimized in DFT corrected by the van der 



Waals interaction as explained 


in the text 










Bond 


a 


h 


c 


d 


e 


0-C 


1.427 


1.428 


1.429 


1.427 


1.427 


C-C 


1.506 


1.509 


1.512 


1.506 


1.506 


C-H 


0.974 


1.107 


1.108 


1.105 


1.106 


0-H 


0.822 


0.982 


0.982 


0.982 


0.982 


Hydrogen 0-H 


1.927 


1.880 


1.868 


1.730 


1.730 


Methyl C-H 


0.987 


1.098 


1.099 


1.098 


1.099 




Figure 1: The labeling scheme for C, H, and O atoms used in Tables 2-5. 
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Table 2: The NMR shieldings a (ppm) in the ethanol crystal, a: atomic positions from X- 
ray measurements all positions are optimized in DFT, c: all positions are optimized 

in DFT corrected by the van der Waals interaction, d: positions of hydrogen atoms are 
optimized in DFT, e: positions of hydrogen atoms are optimized in DFT corrected by the 
van der Waals interaction. Atom numbers are shown in Fig. 1. 



Nucleus 


o(a) 




a(c) 


a{d) 


a(e) 


CI (CH2 group) 


130.2 


109.2 


108.8 


110.7 


110.2 


CIO (CH2 group) 


135.8 


111.9 


111.5 


111.7 


111.2 


C2 (methyl group) 


175.6 


153.9 


152.6 


155.1 


154.6 


Cll (methyl group) 


188.0 


156.9 


155.8 


158.1 


157.6 


H3 (OH group) 


35.6 


30.8 


30.7 


31.4 


31.0 


H12 (OH group) 


36.3 


26.6 


26.3 


25.7 


25.4 


H6 (methyl group) 


31.4 


29.5 


29.2 


30.4 


29.9 


H7 (methyl group) 


32.6 


29.2 


28.9 


29.7 


29.4 


H8 (methyl group) 


33.7 


29.6 


29.3 


30.0 


29.7 


H15 (methyl group) 


34.9 


30.0 


29.7 


30.7 


30.3 


H16 (methyl group) 


33.3 


30.0 


29.8 


30.7 


30.4 


HI 7 (methyl group) 


34.7 


29.7 


29.4 


29.6 


29.1 


H4 (CH2 group) 


31.4 


27.1 


26.9 


27.7 


27.4 


H5 (CH2 group) 


28.6 


26.3 


26.0 


26.8 


26.5 


H13 (CH2 group) 


30.4 


26.5 


26.3 


27.3 


26.9 


H14 (CH2 group) 


30.4 


27.0 


26.8 


27.6 


27.2 


09 


314.9 


255.3 


254.5 


253.0 


252.7 


018 


335.3 


274.5 


273.5 


268.6 


267.7 
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Table 3: The NMR chemical shifts 5 (ppm) in the ethanol crystal, a: X-ray atomic 
positions [1], b: all atomic positions are optimized in DFT, c: all atomic positions are 
optimized in DFT corrected by the van der Waals interaction, d: positions of hydrogen 
atoms are optimized in DFT, e: positions of hydrogen atoms are optimized in DFT cor- 
rected by the van der Waals interaction. The shieldings in TMS and ethanol molecules 
(calculated here under similar conditions as for the crystal) are taken as reference values. 
Atom numbers are shown in Fig. 1. 



Nucleus 


5(a) 


m 


5(c) 


m 


5(e) 


CI (CH2 group) 


48.0 


69.0 


69.4 


67.5 


68.0 


CIO (CH2 group) 


42.4 


66.3 


66.7 


66.5 


67.0 


C2 (methyl group) 


2.6 


24.3 


25.6 


23.1 


23.6 


Cll (methyl group) 


-9.8 


21.3 


22.4 


20.1 


20.6 


H3 (OH group) 


-5.2 


-0.4 


-0.3 


-1.0 


-0.6 


H12 (OH group) 


-6.0 


3.8 


4.0 


4.7 


5.0 


H6 (methyl group) 


-1.0 


0.9 


1.2 


0.0 


0.5 


H7 (methyl group) 


-2.2 


1.2 


1.5 


0.7 


1.0 


H8 (methyl group) 


-3.3 


0.8 


1.1 


0.4 


0.7 


H15 (methyl group) 


-4.5 


0.4 


0.7 


-0.3 


0.1 


H16 (methyl group) 


-2.9 


0.4 


0.6 


-0.3 


0.0 


H17 (methyl group) 


-4.3 


0.7 


1.0 


0.8 


1.3 


H4 (CH2 group) 


-1.0 


3.3 


3.5 


2.7 


3.0 


H5 (CH2 group) 


1.8 


4.1 


4.4 


3.6 


3.9 


H13 (CH2 group) 


0.0 


3.9 


4.1 


3.1 


3.5 


H14 (CH2 group) 


0.0 


3.4 


3.6 


2.8 


3.2 


09 


-56.3 


3.3 


4.1 


5.6 


5.9 


018 


-76.7 


-15.9 


-14.9 


-10.0 


-9.1 
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Table 4: Quadrupolar coupling constants Cq (in MHz) of the ethanol crystal, a: X- 
ray atomic positions b: all atomic positions are optimized in DFT, c: all atomic 
positions are optimized in DFT corrected by the van der Waals interaction, d: positions 
of hydrogen atoms are optimized in DFT, e: positions of hydrogen atoms are optimized 
in DFT corrected by the van der Waals interaction. Atom numbers are shown in Fig. 1. 



l\Fn plpn 


Cr,(o) 


Cr,(h) 








rlo [KJri group J 


U.OU 


U.zo 




n OR 


n OR 

U.ZO 


H12 (OH group) 


0.89 


0.22 


0.22 


0.22 


0.21 


H6 (methyl group) 


0.22 


0.18 


0.18 


0.18 


0.18 


H7 (methyl group) 


0.30 


0.17 


0.16 


0.17 


0.17 


H8 (methyl group) 


0.39 


0.18 


0.18 


0.18 


0.18 


H15 (methyl group) 


0.42 


0.18 


0.18 


0.18 


0.18 


H16 (methyl group) 


0.29 


0.18 


0.18 


0.18 


0.18 


H17 (methyl group) 


0.44 


0.19 


0.18 


0.18 


0.18 


H4 (CH2 group) 


0.46 


0.18 


0.17 


0.18 


0.17 


H5 (CH2 group) 


0.25 


0.15 


0.15 


0.16 


0.15 


H13 (CH2 group) 


0.38 


0.16 


0.16 


0.17 


0.17 


H14 (CH2 group) 


0.36 


0.17 


0.17 


0.17 


0.17 


09 


10.92 


11.47 


11.44 


11.50 


11.48 


018 


10.47 


10.92 


10.90 


10.56 


10.57 



Quadrupolar parameters are shown in Table 4 and Table 5. Although 
existing NQR experiments do not provide signs of quadrupolar coupling con- 
stants we, in principle, show them for calculated quadrupolar parameters. 
However, all calculated quadrupolar parameters are positive. Coupling con- 
stants for H decrease, and for O atoms increase in DFT optimized models. 
Results for all coupling constants and asymmetry parameters in the DFT- 
optimized models b — e are the same or similar to each other. The largest 
differences between the model a and DFT-optimized models b — e exist for 
coupling constants of hydrogen atoms in the OH group. 

Infrared spectra for crystal structures a, b and d are shown in Fig. 2. We 
have calculated infrared spectra for all five models. However, the results show 
that small changes in coordinates induced by the van der Waals corrections 
(models c and e) do not produce substantial changes in infrared spectra. The 
model a shows many (~ 28%) negative vibrational frequencies. In contrast 
to the models b, c, d, and e, the number of negative frequencies in the model 
a does not decrease with changes of parameters typically used in studies of 
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Table 5: Quadrupolar asymmetry parameters rjQ of the ethanol crystal, a: X-ray atomic 
positions b: all atomic positions are optimized in DFT, c: all atomic positions are opti- 
mized in DFT corrected by the van der Waals interaction, d: positions of hydrogen atoms 
are optimized in DFT, e: positions of hydrogen atoms are optimized in DFT corrected by 
the van der Waals interaction. Atom numbers are shown in Fig. 1. 



Nucleus 




VQib) 




VQ{d) 




H3 (OH group) 


0.10 


0.14 


0.14 


0.15 


0.14 


H12 (OH group) 


0.08 


0.16 


0.16 


0.16 


0.16 


H6 (methyl group) 


0.04 


0.03 


0.03 


0.04 


0.04 


H7 (methyl group) 


0.02 


0.02 


0.02 


0.03 


0.03 


H8 (methyl group) 


0.02 


0.02 


0.02 


0.03 


0.03 


H15 (methyl group) 


0.02 


0.02 


0.02 


0.03 


0.03 


H16 (methyl group) 


0.03 


0.02 


0.02 


0.03 


0.03 


HIT (methyl group) 


0.02 


0.03 


0.03 


0.03 


0.03 


H4 (CH2 group) 


0.05 


0.08 


0.08 


0.08 


0.08 


H5 (CH2 group) 


0.06 


0.07 


0.07 


0.08 


0.08 


H13 (CH2 group) 


0.05 


0.06 


0.06 


0.05 


0.05 


H14 (CH2 group) 


0.05 


0.06 


0.06 


0.05 


0.05 


09 


0.75 


0.81 


0.81 


0.82 


0.81 


018 


0.98 


0.88 


0.88 


0.92 


0.91 
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(cm ■'■) 

3600. 3200. 2800. 2400. 2000. 1600. 1200. 800.0 400.0 0.000 




V<fevenumber [cm'-'-) 

3600. 3200. 2800. 2400. 2000. 1600. 1200. 800.0 400.0 0.000 




V\fevenurnber (cm"-'-) 

3600 3200 2800 2400 2000 1600 1200 800, 400. 0.00 




Figure 2: Infrared spectra of ethanol crystal calculated by density functional perturbation 
theory methods, (a) model a: X-ray atomic positions (b) model b: all atomic positions 
are optimized in DFT, (c) model d: positions of hydrogen atoms are optimized in DPT. 
Model a is unstable (see text). 
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vibrational properties of materials by density functional perturbation theory 



methods [26|, |27| . These parameters are accurate self-consistency thresholds 
in the ground state and phonon calculations, as well as the larger energy 
cutoffs. The persistence of many negative frequencies only in the model a, 
indicates its inherent instability. Therefore, as for the bond lengths, the cal- 
culation of vibrational properties of five ethanol crystal models shows that 
new measurements of atomic positions in the ethanol crystal are necessary. 
Spectra in Fig. 2 (b) and Fig. 2 (c) show lines in two intervals of wavenum- 
bers: (500-1500) cm~^ and (2800-3500) cm~^, where IR peaks have also been 
measured jii], 23||. In contrast, the spectrum of model a (Fig. 2 (a)) does 
not have lines in the interval (2800-3300) cm~^. Several far infrared lines in 
Fig. 2 (b) and Fig. 2 (c) also exist, in good agreement with measured data: 



(78-279) cm"^ |22 



4. Conclusions 

Using the GIPAW and PAW DFT methods we calculate NMR chemical 
shielding tensors and NQR parameters of the ethanol crystal. To the best 
of our knowledge, corresponding experimental NMR and NQR data are not 
available. We consider several sets of atomic positions: experimental from 
X-ray diffraction measurements, optimized in the GGA DFT model, and op- 
timized in the GGA DFT model corrected by the van der Waals interaction. 
Situations where only hydrogen atoms are allowed to move, or where all 
atoms move, are studied in the optimization procedure. Differences in atom 
positions between DFT-optimized models and rather old X ray measurements 
are found. It is also calculated that shieldings substantially decrease when 
all atom positions are relaxed in density functional models. Results for in- 
frared spectra of crystal models where positions of either hydrogen, or all 
atoms are optimized, show good agreement with measured data. In contrast, 
the crystal model where all coordinates are taken from X ray measurements, 
in our calculations shows instability, as well as the absence of intervals of 
measured IR lines. New measurements of atomic positions in the ethanol 
crystal are desirable. We hope that our results will be helpful for analysis 
and assignments in future measurements of NMR and NQR parameters for 
the ethanol crystal. 
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